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1. Introduction 



Lattice QCD has played an important role in study of the non-perturbative aspects of QCD. 
However, its application to the finite density system has not been established due to serious diffi- 
culty of the sign problem. In this report we propose a new framework of investigating the 2-flavor 
QCD with finite temperature and density by using the Karsten-Wilczek (KW) lattice fermion ||l|], 
which possesses only two species doublers, i.e. minimally doubled fermion. This lattice formula- 
tion lifts degeneracy of 16 species without breaking its chiral symmetry by introducing a species- 
dependent imaginary chemical potential, instead of a species-dependent mass term introduced in 
the Wilson fermion formalism. Because of the chemical potential term, its discrete symmetry is 
not sufficient to be applied to fully Lorentz symmetric system, i.e. zero temperature and density, 
but enough to study the in-medium QCD. To show the usefulness of the KW fermion, we study 
strong-coupling lattice QCD with temperature and density. 



2. Symmetry of KW-type minimally doubled fermion 

The KW fermion is a kind of minimally doubled fermions, involving only two species dou- 
blers by introducing a species-dependent imaginary chemical potential, which we call "flavored 
chemical potential". This term preserves its chiral symmetry and ultra-locality but breaks 
some of discrete spacetime symmetries [^, 0]. We then need three counter terms to take a cor- 
rect Lorentz symmetric continuum limit: dimension-3, YiJaY = iW^W^ ^^'^ dimension-4 terms, 
\j/Y4d4Y^ FjaFja. [^]. The fermionic part of KW fermion action with the counter terms is given by 



1 4 _ r ^ _ / 



(2.1) 



The second term in the first line including i'Y4 is the flavored chemical potential term, which we also 
call Karsten-Wilczek(KW) term. We here introduce a parameter r in analogy to Wilson fermion. 
/i3 and d4 are parameters for the dimension-3 and dimension-4 counter terms, respectively. The 
corresponding Dirac operator in the momentum space yields 

4 3 

aDKwip) = ' 52 7a( sinap^ + irY4 51 (1 - cosap^) + iHsjA + ''3^474 sinap4, (2.2) 



which has only two zeros at p = (0,0,0, ^arcsin ) when —I — d4 < jJ.^ < \ + d4 with 

r = I. When we expand the Dirac operator around the zeros, its dispersion relation is not Lorentz 
symmetric. As shown in Ref.|Q], the tuning condition for the correct dispersion relation is given 
by {\ +d4)^ = 1 + jLt| at the tree level. Moreover, it is shown that /I3 has to be tuned to control 
imaginary chemical potential in ff{l/a). 

Symmetries of the lattice action ( p] ) are chiral symmetry, cubic symmetry corresponding to 
permutation of spatial three axes, CT and P [||] : 
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1. U (1) chiral symmetry (75 ® T3 [|, 0, §]) 

2. Cubic symmetry 

3. CT 

4. P 



It is notable that these symmetries are the same as those of the finite-density lattice QCD: As an 
example, we look into the naive lattice action with chemical potential, which is given by 



1 



L ^xjj (u,^x+jVx+J- Ux,x-jVx-j) + VxY4 (e^U^,,+iVx+4 - ^ ^U^,x-4Vx-4 
7=1 



(2.3) 

The 4th direction hopping term, involving chemical potential, breaks the hypercubic symmetry into 
the spatial cubic symmetry, and also C, P, and T into CT and P, which are the same symmetries of 
(p|). It means that, even if we introduce chemical potential as (2.3) to KW fermion, the symmetries 
are unchanged. The KW fermion with the exponential form chemical potential is given by. 



9 L VxYj ( U,,,+jVx+} - U,,,^j¥x-}) + h L VxYa [2Vx - U,^x+j¥x+} " ^x.x-j^x-] 
^ j=i ^ ^ j=i 



l+d4 



¥xYA ( U^,^^^ Yx+4 ^ Ux.x-4 Vx-4 ) + ¥xY4 Vx 



(2.4) 



From the viewpoint of the universality class, these two theories, ( |2.3[ ) and ( |2.4| ), should belong to 
the same class. 

Here we remark the way of introducing chemical potential. It was pointed out in [Q] that a 
naive form of the chemical potential, ^y^W = I^VY^V' violates the Abelian gauge invariance and 
requires a counter term to make thermodynamical quantities finite. On the other hand, in the KW 
fermion, the flavored chemical potential term is introduced in this naive form. It leads to necessity 
of tuning to deal with if?" {I /a) additive renormalization of chemical potential, as with the mass 
renormalization in the Wilson fermion. This renormalization effect is relevant to the phase diagram 
in the (jU3-g^) parameter plane, as discussed in 



3. Strong-coupling lattice QCD 



We study QCD phase diagram in the framework of the strong-coupling lattice QCD with 
KW-type minimally doubled fermion. We extend the strong-coupling analysis with this lattice 



fermion [ ]10[ ] to the finite temperature and density system [ ]11| , |12| , |13[ ]. The effective potential in 
terms of the meson field is obtained by performing the 1-link integral in the strong coupling limit 
(g2 _^ ^-j^ ^j^j jj^gj^ introducing auxiliary fields to eliminate the 4-point interactions. In the case 
with KW fermion, we have to consider both of the scalar a = {yj/Y) ^rid vector 714 = ( 1/^/74 
condensates. Identifying 1/^74 = y"^' the latter corresponds to the imaginary density {(y'^y)- For 
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the case with SU {Nc) gauge group and d = D+l dimensions in the finite temperature and density, 
we obtain the the following effective potential, 

^ss{a,7i4;m,T,i^,l^3,d4) = ^ {{l + r^)o^ + {l-r^)nl) -A^^logA- ^log I ^ det(2„+,_^-)i<i,y<A'. ) 

(3.1) 

In particular, the determinant part for A^^. = 3 is given by 

£det(e«+,--y)i<,-j<;V^ 

8( l + 12cosh2- + 8cosh'^- ) ( 15 - eOcosh^ - + 160cosh^ - - 32cosh^ - + 64cosh*^ - 



-64cosh^cosh- ( -15 +40cosh-' - + 96cosh'* - + 320cosh^ - 



2.11 f K E E 

+ 80 cosh — 1 + 6 cosh^ - + 24 cosh"^ - + 80 cosh*" - 
r V T T T 

+80cosh^cosh| (^-l+cosh^^j +2cosh^, (3.2) 
with 

£ = arcsinh(B/A), = {l+d4f + (^^3 + Dr -^{l - r^)K4^ , B = m+^{l+t^)o. (3.3) 

Here the baryon chemical potential is defined as plb = ^pi- Remark that the next-leading order 
terms in 0'{\/^/D) are omitted in the derivation. See [ ]T4| ] for the detailed calculation. 

In the zero temperature case, we can solve the equilibrium condition analytically. For D = 3 
{d = 4) with m = and r = 1 the potential is given by 

=^eff(a) = ^a^-^log ((1 + ^4)2 + (^3 + 3)2) - max J 3 arcsinh ( -—=^==== 

2 ' y \^^(i+cf4)2+(^3+3)2y J 

(3.4) 

In this case there are two local minima of the free energy as a function of a at a = and a = Oq. 
This Oq can be determined by the gap equation, d^^s/do 



= 0, 

CT=f70 



(1+^4)' + (M3 + 3)' 

18 



((l+^/4)2 + (M3 + 3)2)2 



(3.5) 



Comparing these two local minima, we can show that the global minimum changes from o = Oq 
to a = at the critical chemical potential as 

critical ( J = 0) = 3 arcsinh ( 1 - -a,?- (3.6) 

This chiral phase transition is of 1st order because the order parameter a changes discontinuously 
at this critical chemical potential. We can also evaluate the baryon density ps = — 5=^eff/^MB at 
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Figure 1: (Left) Phase diagram for the chiral transition with r — \, pL^ = —0.9 and cIa, = 0. Green and red 
Unes show 2nd and 1st transition lines, respectively. The transition order is changed from 2nd to 1st at the 
tricritical point {jx^ ,T^) = (0.804,0.234). (Right) Three-dimensional chiral phase diagram for T , pig and 
/i3 for m = where /is runs within half of the physical range — 3 < /X3 < ^IfLjl — 3. Green, red and purple 
lines show 2nd, 1st order transitions and tricritical point, respectively. 



r = 0. It turns out to be empty Pe = when < ;U|"°'''''. On the other hand, when \Ib > jUg 
it is saturated as = 1. 

We then discuss the phase diagram with respect to chiral symmetry. We now concentrate on 
the case with r = 1 for simplicity because the effective potential (3.1) is independent of 714 in such 
a case. The 2nd order chiral phase boundary is given by the condition, such that the coefficient of 
in the effective potential ( 3J_ ) becomes zero. When the order of the phase transition is changed 
from 2nd to 1st, the coefficient of as well as should vanish. The left panel of Fig. |] shows the 
phase boundary of the chiral transition with r = 1, ^1^3 = —0.9 and m = for = 0. The counter 
term parameter is taken from the physical region -73277 <Ai3 + 3 < 73277 [Hi. The order of the 
phase transition is changed from 2nd to 1st at the tricritical point (/Xg',r'") = (0.804,0.234). We 
also depict a condensate and the baryon density pe = —d^^^/dpLB as functions of }JLb with several 
fixed T in Fig. ^ We find that there are 1st {T < T'") and 2nd {T > T^^) order phase transitions for 
a, followed by the phase transition of the density ps- For m / 0, we can show that the crossover 
transition instead appears with the 2nd order critical point. 

These results are qualitatively consistent with those with strong-coupling lattice QCD with 
staggered fermions, while there are some quantitative differences. For example, the KW phase 
diagram is suppressed in T direction compared to that in staggered. We here compare the ratio 
of the transition baryon chemical potential at T = to the critical temperature at /ig = 0, = 
H^(T = 0)/Tc{hb = 0). In staggered fermion, this ratio is /?° - 3 x 0.56/(5/3) ~ 1 [|l|, |l3]], while 
R% ~ 0.767/0.356 ~ 2.2. In the real world, this ratio is larger, > Mn/HO MeV ~ 5.5. When 
the finite coupling and Polyakov loop effects are taken into account for staggered fermion, Tc{ij.b = 
0) decreases, = 0) stays almost constant, then /?" value increases [|r|]. Larger R° with KW 
fermion in the strong coupling limit may suggest smaller finite coupling coiTcctions in the phase 
boundary. Another interesting point is the location of the tricritical point. In KW fermion, the ratio 
is /Jjriy = 0.804/0.234 ~ 3.4, while 7?^ = 1.73/0.866 ~ 2.0 for unrooted staggered fermion [[12, 
13]. It would be too brave to discuss this value, but R^£^ is consistent with the recent Monte-Carlo 
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m = 0, r = 0.3 771 = 0, r = 0.2 




Baryon Chemical Potential Baryon Chemical Potential 



m = 0.1, T = 0.3 m = 0.1, T = 0.2 




Baryon Chemical Potential Baryon Chemical Potential 



Figure 2: Chiral condensate O and the baryon density ps for (left) T = 0.3 and (right) T ~ 0.2 with di^ ~ 0. 
Top and bottom panels show the massless m = and massive m = 0.1 cases. There are 1st and 2nd phase 
transitions for a. In the case of m ^ 0, there appears the crossover behavior instead of the 2nd order 
transition. 



simulations (see references in [[T^), which impUes that the critical point does not exist in the low 
baryon chemical potential region, jUg/r < 3. These observations reveal usefulness of KW fermion 
for research on QCD phase diagram. 

Apart from the phase transitions, the jig dependence of o and ps seems to have some charac- 
teristics in Fig. ^ At r = 0.3 > T'" with m=0, a and undergoes the 2nd-order phase transition 
at jiB ^ 0.5, and at a larger jUg (/ig ~ 1.15), increasing rate of pe as a function of /ig becomes 
higher again. At lower temperature, T = 0.2 < T'", partial restoration of the chiral symmetry is 
seen before the first order phase transition. Since we have not taken care of the diquark condensate, 
these continuous changes are not related to the color superconductor. Other types of matter, such 
as quarkyonic matter [17], partial chiral restored matter [jlsl], or nuclear matter, may be related to 
the above chai^acteristics. 

In this report we focus on the case with r = I and = for simplicity. It is shown in [ |l4| ] that 
effects of these parameters are just quantitative. 



4. Summary 

We have proposed a new framework for investigating the two-flavor finite-(r,/x) QCD phase 
diagram. We have shown that the discrete symmetries of KW fermion strongly suggest its appli- 
cability to the in-medium lattice QCD. To support our idea, we study the strong-coupling lattice 
QCD in the medium and derive the phase diagram of chiral symmetry for finite temperature and 
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chemical potential. We have obtained the phase diagram with 1st, 2nd-order and crossover critical 
lines, which is qualitatively in agreement to results from the model study. 
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